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Abstract 

We have measured scintillation properties of pure Csl crystals used in the shower 
calorimeter built for a precise determination of the 7r + — >7r°e + f e decay rate at the 
Paul Scherrer Institute (PSI). All 240 individual crystals painted with a special 
wavelength-shifting solution were examined in a custom-build detection apparatus 
(RASTA — radioactive source tomography apparatus) that uses a 13 Cs radioactive 
gamma source, cosmic muons and a light emitting diode as complementary probes 
of the scintillator light response. We have extracted the total light output, axial 
light collection nonuniformities and timing responses of the individual Csl crystals. 
These results predict improved performance of the 3ir sr PIBETA calorimeter due 
to the painted lateral surfaces of 240 Csl crystals. The wavelength-shifting paint 
treatment did not affect appreciably the total light output and timing resolution of 
our crystal sample. The predicted energy resolution for positrons and photons in the 
energy range of 10-100 MeV was nevertheless improved due to the more favorable 
axial light collection probability variation. We have compared simulated calorimeter 
ADC spectra due to 70 MeV positrons and photons with a Monte Carlo calculation 
of an ideal detector light response. 
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1 Introduction 



The PIBETA detector at the Paul Scherrer Institute in Switzerland was de- 
signed for the precise measurement of the pion beta {it (3) decay branching ratio 
(7r + ^7r°e + z/ e ) as well as for study of other rare pion and muon decays [1]. The 
experiment has started to take production data in the summer of 1999. The 
heart of the detector is the ~37rsr pure cesium iodide (Csl) segmented calo- 
rimeter. In-depth knowledge of the calorimeter response to the positrons and 
photons in the energy range 10-100 MeV is essential in order to achieve the 
proposed accuracy of the branching ratio measurements of ~0.5 %. Experimen- 
tal signature of a tt/3 event is determined by the prompt decay tt° — > 77. The 
(1.025± 0.034) xlO -8 pion beta decay branching ratio [2] will be remeasured 
relative to the 10 4 times more probable ir + ^e + v e decay rate, that is known 
with the combined statistical and systematic uncertainty of ~0.33% [3,4]. 

The detector energy response functions were determined both experimentally, 
by scanning the calorimeter with the monoenergetic e + and 7 beams, and in 
Monte Carlo simulations. Inputs to the simulation calculations were the aver- 
age light outputs and light output nonuniformity functions for each individual 
calorimeter module. In this paper we describe how these input parameters, 
essential for a realistic Monte Carlo detector simulation, have been measured 
in a dedicated tomography apparatus. Section 2 details the surface treatment 
of the Csl scintillator crystals. Design and operation of the radioactive source 
scanning apparatus are explained in Sec. 3. The extracted light output param- 
eters are discussed in Sees. 4 and 5. In Sec. 6 we explain how these parameters 
are used in a GEANT-based calculation of the calorimeter response. 



2 Csl Crystal Surface Treatment 

The PIBETA calorimeter consists of 240 truncated crystal pyramids with a 
height of 22 cm, which corresponds to 12 radiation lengths. Nine different 
detector shapes were required to complete the crystal sphere [5]: four ir- 
regular hexagonal truncated pyramids (we labeled them HEX-A, HEX-B, 
HEX-C, and HEX-D), one regular pentagonal (PENT) and two irregular 
half-hexagonal truncated pyramids (HEX-D 1 and HEX-D2), and two trape- 
zohedrons (VET-1, VET-2). The volumes of studied Csl crystals vary from 
797cm 3 (HEX-D1/2) to 1718cm 3 (HEX-C). 

Twentyfive pure Csl crystals were produced by Bicron Corporation facility in 
Newbury, Ohio. The reminder of the crystals were grown and cut in the In- 
stitute for Single Crystals (AMCRYS) in Harkov, Ukraine. The more detailed 
description of the PIBETA calorimeter design is given in Refs. [6] and [7]. 
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In Fig. 1 we show the partially assembled calorimeter with front faces and 
wrapped side surfaces of exposed Csl modules visible inside the mechanical 
support structure. 

The light output of each crystal, namely the average number of photoelectrons 
for unit energy deposition caused by ionizing charged particles, can be tuned 
by using the different polishing, matting and wrapping treatments of the de- 
tector surfaces. The goal of such treatments is usually to maximize the light 
yield thus obtaining the best energy resolution, while keeping the detector 
response linear irrespective of the location of the energy deposition. 

Our initial treatment of the crystal surfaces involved a test sample of 74 Csl 
detectors and consisted of: (a) polishing the crystal surfaces with a mixture 
of 0.2 jum aluminum oxide powder and etylenglycol, (b) wrapping the lateral 
surfaces with two layers of Teflon foil plus one layer of aluminized Mylar sheet, 
and (c) covering the front crystal surface with a black paper template. The 
crystal back surface with the glued phototube was left uncovered. We reported 
the predicted response of the calorimeter with all Csl modules treated in the 
described manner in our previous paper [7]. 

Further test measurements and the technical developments [8] suggested an 
improved optical treatment of crystal surfaces. Due to the fact that Csl is 
slightly hygroscopic, it is desirable to "seal" the scintillators against poten- 
tial surface deterioration and thus minimize changes of the detector response 
throughout the duration of the experiment. Traditional polishing and wrap- 
ping methods cannot ensure long term stability of the relevant scintillation 
parameters due to the insufficient chemical resistance of alkali halides [9]. 
Furthermore, a suitable protective coating of the crystal surfaces can at the 
same time be used to modulate the scintillator light collection properties and 
can be easily removed, if necessary. 

We therefore decided to paint the surfaces of every Csl crystal using a spe- 
cial organosilicon mixture. The solution was developed at the Harkov Sin- 
gle Crystals Research Institute [10]. For our purposes the most useful poly- 
mers were polyphenylsiloxanes due to their high thermal stability, good ra- 
diation resistance and the moisture-proof quality. The optical properties of 
the organosilicon coatings depend on the nature of the organic radicals. Their 
phenyl polymer bonds can efficiently transfer the energy to the luminescent 
admixtures [11]. The selected polymers and solvents have a nonpolar character 
resulting in a low moisture absorption. On the other hand, the high optical 
transparency of the varnishes was achieved in a purification process by the 
absorption trapping of silica, aluminum oxide and hydroxide impurities. The 
resulting absorption and reemission of light by the coating in the wavelength 
range above 300 nm was better than 99%. 
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Specifically, the Waveshift lacquer was the ladder organosilicon copolymer with 
the chemical composition PPO+POPOP+COUM.l, where PPO is 2.5-Di- 
phenyloxazole, POPOP represents 1.4-Di-2-(5-Phenyloxazolile-Benzene) and 
COUM.l is 7-Diethylamine-4-Methylcoumarin [12]. That paint solution was 
applied twice resulting in a protective cover with a total thickness of 120 /zm. 
The dried paint layer has the index of refraction tid=1.62 and the density 
p=1.094, and is soluble in ethyl alcohol. Following the paint application each 
detector module was wrapped in two layers of a 38 /im PTFE Teflon sheet 
and two layers of 20 /im thick aluminized Mylar to make it optically insulated 
from its calorimeter neighbors. 

We have compared the average fast-to-total (F/T) analog signal ratio and 
the photoelectron statistics of all lacquer-treated crystals with the measure- 
ments of polished and Teflon- wrapped crystals we reported on in Ref . [7] . The 
F/T ratio was virtually unchanged at 0.77, while the number of photoelec- 
trons per MeV was just under 3% lower, changing from an average of 64 to 
62 Photoelectrons/MeV at the ambient temperature of 18°C. The combined 
statistical and systematic uncertainty of these measurements was less then 1 
Photoelectron/MeV. 

The averaged lineshape of an analog cosmic muon signal from one Csl detector 
captured on a digital oscilloscope is shown in Fig. 2. The 10-to-90 % amplitude 
rise time of the analog signal changed from 4.0 ns for the Teflon- wrapped 
scintillators to 4.8 ns for the lacquer-painted detectors, due to the fact that 
a fraction of the scintillator light is now absorbed by the wavelength-shifting 
surfaces and reemitted later with the different wavelength. Histograms of the 
F/T ratios and the photoelectron numbers for the painted crystals are shown 
in Figs. 3 and 4, respectively. 



3 RASTA: RAdioactive Source Tomography Apparatus 

The total light output and the position-dependent light collection probabilities 
for each Csl detector were measured in a test setup using a 137 Cs radioactive 
gamma source with an activity of 3.7 MBq. It was found that this method 
produced results consistent with the cosmic muon tomography [7] and was well 
suited for the routine crystal uniformity tests. Thus an apparatus was designed 
and built to allow a fast, precise and reproducible method of determining the 
optical nonuniformity of each Csl crystal. 

The drawing of the main parts of the apparatus is shown in Fig. 5. The 137 Cs 
radioactive gamma source was embedded in a lead collimator. The collima- 
tor had a thickness of 5 cm while its pin-hole had a diameter of 6 mm. The 
probability for the 0.662 MeV 137 Cs photons to penetrate that lead thickness 
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is < 1 %. The ionization energy deposited at the central axis of a Csl crys- 
tal is about one-sixth of the energy deposition near the crystal surface. The 
collimator was mounted on an aluminum plate, which was moved by a step- 
ping motor. The precision of the step motor movements was 12.5 /im. The 
examined crystal was placed on a horizontal aluminum plate using a dedi- 
cated positioning system with a precision of ~0.3mm. The stepping motor 
was controlled by a PC486 personal computer through a RS232 interface. 

The schematic diagram of the logic electronics is shown in Fig. 6. The analog 
photomultiplier signal was amplified by a factor of 10 using an LRS 612A 
photomultiplier amplifier. The signal was then split, one copy going into the 
trigger logic while the other copy was delayed and fed into an LRS 2249A 
CAMAC ADC. The trigger created a computer LAM in the SIN IO506 in- 
put/output unit. The CAMAC modules were controlled using a HYTEC 1331 
Turbo CAMAC-PC interface. 

The stability of the system after switching on the PMT high voltage was 
studied in preliminary calibrations. The data were taken with a Csl crystal 
in a fixed position every 5 minutes and the results are shown in Fig. 7. We 
concluded that the PMT high voltage should be switched on for at least one 
hour before doing any measurements with the detector. 

A full calibration was performed on each Csl crystal by scanning it four times 
within three hours. The measured light output as a function of scan posi- 
tion are shown in Fig. 8 for the Csl detector number S040. After the cross- 
normalization the results of different measurements agreed within the 0.5 % 
statistical uncertainty. 



4 Total Light Output of the Csl Detectors 

Scintillation light of a full-sized Csl crystal was read by a 78 mm diameter 
EMI 9821QKB photomultiplier tube with quartz photocathode window. The 
smaller half-sized crystals were directly coupled to a 52 mm EMI 9211QKA 
PMTs with similar characteristics [13]. Both phototube types were glued to 
the back faces of the crystals using a clear silicone elastomer Sylgard 184 
manufactured by Dow Corning [14]. 

A light-emitting diode (LED) was attached to the exposed back side of the Csl 
detector. The whole detector assembly was placed in the light-tight RASTA 
box. The LEDs were excited at a ~10 Hz rate using a NIM diode driver with an 
adjustable output voltage. The LED pulses simulated constant energy depo- 
sitions in the Csl crystal corresponding to a 10-100 MeV range. The absolute 
energy calibration was established by comparison with measured cosmic muon 
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ADC spectra. These spectra were acquired in calibration runs triggered by two 
small plastic scintillator tag counters placed above and below the Csl crystal. 
Their expected peak positions, expressed in the units of MeV, were calculated 
in a GEANT Monte Carlo simulation. 

For LED measurements the variance a\ of a PMT pulse-height spectrum 
depends upon the mean number of photoelectrons N pe from the photocathode 
created per unit energy deposition in the scintillator: 

o% = T,°i+E/*p» (!) 



where E is the spectrum peak position and of 's are assorted variances, such as 
the instabilities of LED driving voltage and temporal ADC pedestal variations. 

The energy resolutions erg's were measured for 9 different values of the LED 
driving voltage and the photoelectron numbers N pe were calculated by the 
least-squares fit using Eq. (1). The mean numbers of photoelectrons normal- 
ized to the unit equivalent-energy deposition in MeV and scaled down the 
ambient temperature of 18°C are shown in Fig. 4. The histograms for the full- 
sized crystals and half-sized crystals are displayed separately, with the ratio 
of 2.4 in photoelectron statistics essentially consistent with the surface area 
ratio of 2.25 in the two differently-sized photocathodes. 



5 Spatial Light Output Nonuniformity of the Csl Detectors 

The 3-dimensional spatial distribution of scintillation light output can be spec- 
ified by giving the number of photoelectrons N pe (x, y, z) produced by a 1 MeV 
energy deposition at the point (x,y,z) ("3-dimensional light nonuniformity 
function"). In the following discussion we limit ourselves to the linear and 
one-dimensional variation of the detected light output N pe (z) along the long 
axis of the detector: 

f N 1 + a zl - z < z < 10 cm, 
Npe(z) oc < N 2 + a z2 -z 10 < z < 18 cm, (2) 
I N 3 + a z3 ■ z 18 < z < 22 cm 

where a 2 's are the linear optical nonuniformity coefficients in the axial coor- 
dinate z, and the coordinate system origin is at the center of the detector 
front face. In this approximation the transverse light output variations are 
integrated away. In our figures and tables we are quoting the a z light nonuni- 
formity parameters in units of %/cm. 
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The complete RASTA database of all 240 studied Csl crystals is made avail- 
able for online inspection at the PIBETA WWW site [15]. The measurements 
for each Csl crystal are summarized on separate WWW pages. Individual 
pages start with the pulse-height spectra recorded with the cosmic muon trig- 
ger that is used for the absolute energy calibration of the ADC scale. The 
LED "energy" spectra corresponding to 9 different LED intensities are shown 
next, together with the superimposed Gaussian fits. The least-square linear 
fits of the data using the Eq. (1), which allow the deduction of the average 
number of photoelectrons per MeV, as well as the associated fit uncertainties, 
are documented below the spectra. The timing resolution histogram is also 
presented; it gives the standard deviation of the Csl detector timing spectrum 
with respect to the tag counters TDC value for the cosmic muon triggers. 

Also displayed are the ADC spectra taken with the 137 Cs radioactive source 
at 11 equidistant positions along the long axis of a studied crystal, together 
with Gaussian-exponential fits. The fitted Gaussian peak positions are plotted 
in the axial light nonuniformity figures that follow next. The relative axial 
variations of the light collection nonuniformities are shown for 6 representative 
Csl crystals in Fig. 9. The extracted linear nonuniformity coefficients a z 's are 
listed last. 

The relation between the linear nonuniformity coefficients a z \ in the front part 
of the crystal (z < 10 cm) and the coefficients a Z 2 describing the position- 
dependent change of the axial light output in the central and back crystal 
part is shown on the scatter plot in Fig. 10 for all 240 crystals. Fig. 11 com- 
pares the linear function a Z 2 = c\ + c 2 • a z \ fitted to the scatter plot with 
the result deduced from our previous measurements of the Teflon- wrapped 
crystals, Ref. [7]. The parameters c\ and c 2 were determined by the "robust" 
estimation method [16] by imposing the requirement of the minimum abso- 
lute deviation between the measured and calculated values. The treatment of 
crystal surfaces with the wavelength-shifting lacquer produces more favorable 
light output collection along the detector axis by increasing the collected light 
from the back part of a crystal and thus compensating for the energy leak- 
age of electromagnetic showers and improving the detector energy resolution. 
These results were confirmed in the full-fledged GEANT simulation described in 
Sec 6. 

In Fig. 12 we show the timing resolutions for the all measured Csl crystals. We 
have histogrammed the standard deviations of the Gaussian functions fitted 
to the spectrum of the timing differences between the Csl detector and the 
small scintillator tag counter just above the crystal. The average Csl detector 
timing uncertainty specified in such a way is 0.68 ns. 

We have summarized the average scintillation properties of all 240 studied Csl 
crystals in Tables 1 and 2. 
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6 Simulated ADC Lineshapes of 70 MeV/c e ±, s and 7's 

We have developed a comprehensive GEANT Monte Carlo description of the 
PIBETA detector that includes all major active as well as passive detector 
components [17]. The user code is written in a modular form in standard 
FORTRAN and organized into over 300 subroutines and data files [18]. The 
graphical user interface to the code based on Tel routines supplements the 
program. 

Standard GEANT routines account for electromagnetic shower processes in Csl: 
leak-through, lateral spreading, backsplash, as well as the electromagnetic in- 
teractions in the active target and tracking detectors: bremsstrahlung, Bhabha 
scattering, in-flight annihilation for positrons, and pair production, Compton 
scattering, etc. for photons. 

The individual PIBETA detector components, both passive and active ones, 
can be positioned or "switched off" without recompiling the code. The follow- 
ing detectors are defined in the Monte Carlo geometry: 

(1) a beam counter, an active plastic scintillator degrader and several versions 
of segmented stopping active targets; 

(2) two concentric cylindrical multiwire proportional chambers used for 
charged particle tracking; 

(3) a 20-piece cylindrical plastic scintillator hodoscope used for charged par- 
ticle detection and discrimination; 

(4) a 240-module pure Csl calorimeter sphere; 

(5) a 5-plate cosmic muon veto scintillator system; 

(6) passive calorimeter support structure and the individual detector support 
systems, 

(7) the detector phototubes and HV divider bases, and the lead brick shield- 
ing structure. 

User input to the simulation code requires the detector version, e.g., the run 
year, the beam properties, the selected reaction/decay final states, and the 
version of the ADC and TDC simulation codes. User can select one partic- 
ular final state or any combination of different final states with the relative 
probabilities defined by the reaction cross sections and decay branching ratios. 
The simulation of ADC values and TDC hits takes into account the individ- 
ual detector photoelectron statistics, the axial and transverse light collection 
detector nonuniformities, ADC pedestal variations, electronics noise, and the 
event pile-up effects. The photoelectron statistics and the individual detec- 
tor light collection nonuniformity coefficients are initialized from the RASTA 
database file. Calculation of accidental coincidences up to the fourth order is 
optional. 
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Selectable options also include the effects of: 

(1) lateral and axial extent and the divergences of the stopping pion beam; 

(2) positron and muon beam contamination; 

(3) photonuclear and electron knockout reactions in Csl material [19]; 

(4) aluminized Mylar wrapping of Csl modules and plastic scintillator veto 
staves; 

(5) cracks between the Csl detector modules; 

(6) temperature coefficients of the individual calorimeter modules; 

(7) gain instability of the calorimeter detector modules; 

(8) electronics discriminator thresholds and ADC gate width. 

Particular attention was paid to the correct accounting of "software gains" 
for the individual Csl modules. In the simulation calculation the values of the 
detector software gains allow the same user control element as the detector 
high voltages in the operation of the physical detector. The software gains are 
determined in an iterative procedure constrained by the Monte Carlo 7r — > e + v 
positron ADC spectra in each Csl detector. Namely, in the real experiment 
the high voltages of the individual Csl detectors are also set by matching the 
positions of the 69.8 MeV positron peaks. 

The output of the calculation provides the selected experimental layout and 
the final states chosen in the calculation, as well as the used production cross 
sections and branching ratios. A number of physical variables that could be of 
interest to a user are histogrammed in 1- and 2-dimensional format. Raw- wise 
and column-wise PAW Ntuples [20] are used to digitize the individual events. 
Simulated energy depositions and the ADC and TDC values associated with 
the active detectors are saved on an event-per-event basis. A single event 
display program can be used to examine in detail the particular event in an 
output stream. 

Simulated calorimeter ADC spectra for 69.8 MeV monoenergetic positrons and 
70.8 MeV photons generated in the stopping target are shown in Fig. 13. The 
figure shows the ADC lineshapes summed over all 220 calorimeter signals 
that make up the low threshold (-Ethr>5 MeV) trigger logic. The top panel 
shows the spectra without the energy cut on the 20 summed veto Csl signals 
that define the beam opening/exit calorimeter surfaces, the bottom panel 
shows the same spectra with the imposed veto energy cut. Fig. 14 displays the 
same variables but now summed only over the calorimeter detector with the 
maximum energy deposition and its nearest neighbors. That restricted ADC 
sum should reduce the contribution of accidental event pile-up to the smearing 
of the energy spectra. 

The e + and 7 peak positions, their widths and tail fractions of the simulated 
ADC lineshapes are summarized in Tables 3 and 4. We see that the calorimeter 
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energy resolution is improved and the energy spectra tail fractions reduced in 
the detector with the paint-treated crystals. The improvement is smaller for 
the positron showers which deposit most of their energy close to the incident 
particle impact point, but is considerable, over 15%, for the 70MeV photons 
which penetrate deeper into the calorimeter. 



7 Conclusion 

We have measured the scintillation properties of the individual pure Csl crys- 
tals used in the PIBETA calorimeter. The crystals were treated with a special 
wavelength-shifting organosilicon paint. The fraction of the fast decay compo- 
nent in the total light output, the photoelectron statistics and the temperature 
dependence of the light output are documented for the each detector. A 137 Cs 
radioactive source scanning apparatus was designed and used to determine 
the axial light collection nonuniformity of all scintillator crystals. An individ- 
ual Csl detector response was described by three linear light nonuniformity 
coefficients. 

These experimentally determined scintillation properties of the calorimeter 
components provide the essential input to a comprehensive GEANT Monte Carlo 
of the PIBETA detector. We have compared the simulated performance of 
the PIBETA calorimeter having the waveshift-painted Csl crystals with an 
idealized calorimeter characterized by the uniform light collection from every 
volume element of the scintillator. The comparison shows that the PIBETA 
calorimeter made up of the painted Csl crystals has better energy resolution, 
while its timing resolution is not compromised. 
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9 Appendix A. RASTA Programs 

The main computer program which operated the apparatus was named 
RASTA. Its main tasks were to: 

(1) control the stepping motor operation; 



10 



(2) control and read out the CAMAC modules; 

(3) perform histogram analysis including pedestal and background subtrac- 
tion, histogram filling and fitting; 

(4) provide graphical data display. 

Upon the execution of the program, the operator was first asked for the se- 
rial number Sxxx of the crystal being tomographed. A new disk subdirectory 
hisxxx was created in which all the relevant histograms of a scan process were 
stored. 

The starting parameters of the program were defined in the file RASTA. INI: 

[Scan] 
Time=20 
Origin=l .0 
Stepsize=2 . 
Numberof Steps=ll 
[Fit] 

PeakFindMin=10. 
PeakFindMax=800 . 
Transition=-20 . 
NsigmaXLow=2 . 
NsigmaXHigh=6.0 

The category Scan controlled the scanning process: 

(1) Time is the scan time per point in seconds; 

(2) Origin is the starting point of the scan (with the origin marked in the 
drawing of the apparatus, Fig. 5); 

(3) Stepsize is the distance between measured points given in centimeters. 

The category Fit executed the fitting process: 

(1) PeakFindMin (PeakFindMax) is the lower (upper) limit of the peak posi- 
tion passed to the fitting program when searching for the maximum ADC 
value and FWHM of the pulse-height spectrum; 

(2) Transition is the value of the Gauss-exponential transition of the fit 
function measured from the peak position; 

(3) NsigmaXLow and NsigmaXHigh define the histogram range to be fitted in 
units of the Gaussian standard deviation from the peak position. 

After having changed or accepted the default run parameters the step motor 
was initialized and the collimator moved to the parking position where the 
pedestal data and background data were taken. 
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The crystal was then scanned automatically at the points specified in the 
RASTA . INI file. Fig. 15 shows the computer display as it appeared during the 
scanning process. 

At the end of the scan all histograms were stored so that the data could be 
easily reviewed in a Microsoft EXCEL spreadsheet (histograms numbers and 
1) or with the in-house developed Histogram Manager utility [21] (histograms 
2 to NSTEP+3). 

Two other programs are supplied together with the RASTA package: 

(1) HVSET, which helps select the high voltage (HV) value for the photomul- 
tiplier attached to a new crystal being prepared for tomography. The 
additional input parameters are the HV starting value, the preset pulse- 
height spectrum peak value and maximum differ ence between the preset 
and fitted peak value; 

(2) MOVE, which allows to position the collimator at the required location. 
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Fig. 1. Picture of the Csl calorimeter during the mechanical assembly. About half 
of individual crystals (of the 240 total) are seen in a self-supporting configuration 
inside the detector cradle. 

Fig. 2. Comparison of the cosmic muon waveforms for a hexagonal Csl detector 
wrapped in aluminized Mylar (dashed line) and for the same crystal painted with 
the wavelength-shifting lacquer (full line). 

Fig. 3. Histogram of the fast-to-total light component ratio for all studied Csl crys- 
tals. Light output of unwrapped painted crystals was measured with a digital oscil- 
loscope at the average ambient temperature of 22° C. 

Fig. 4. Histogram of the number of photoelectrons per MeV of deposited energy for 
full-sized (left panel) and half-sized (right panel) Csl detectors normalized to the 
ambient temperature of 18°C. 

Fig. 5. Schematic drawing of the RASTA experimental apparatus. The depicted 
assembly is enclosed in a light-tight aluminum box with a hinged cover plate. 

Fig. 6. Schematic drawing of the RASTA logic electronics. In routine data acquisi- 
tion mode the trigger is provided by the discriminated single Csl detector signal. 

Fig. 7. Plot of the drift and leveling off of the pulse-height spectrum peak position 
due to the settling down of the photomultiplier during the first hour after the PMT 
high voltage is switched on. 

Fig. 8. Plots of the reproducibility of the axial light nonuniformity measurements 
demonstrated with four sets of measurements of a single Csl crystal taken consec- 
utively. 

Fig. 9. Plots of the relative light output as a function of the radioactive source axial 
position for six different Csl detectors. The detector phototubes are positioned at 
z=22 cm. 

Fig. 10. Light nonuniformity coefficients a z \ and a Z 2 (%/cm), determined separately 
for the front (z<10 cm) and back (z>10 cm) crystal sections, are plotted in a 
scatter-plot for 240 Csl detectors. Full-sized hexagonal and pentagonal detectors 
are represented with full circles, while the open circles indicate the half-hexagonal 
and trapezial detector shapes. 

Fig. 11. Plot of the relationship between the linear light output nonuniformity co- 
efficients a z \ and a Z 2 for the polished and Teflon- wrapped Csl crystals (dotted line) 
and for crystals painted with the Waveshift lacquer (full line). 
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Fig. 12. Cosmic muon timing for one representative Csl crystal with respect to the 
plastic scintillator tag counter (top panel) and the timing resolutions of all studied 
Csl crystals (bottom panel). The resolution is defined as the standard deviation of 
the TDC timing difference between the tag counter and the Csl detector. 

Fig. 13. The predicted PIBETA calorimeter spectra for monoenergetic 69.8 MeV 
positrons and 70.8 MeV photons. All Csl values above the 0.5 MeV TDC threshold 
(excluding the 20 trapezial veto counters) were summed to obtain the deposited en- 
ergy. The GEANT simulation used different numbers of photoelectrons/MeV and lin- 
ear axial light collection nonuniformities for each crystal as extracted in the RASTA 
analysis. The bottom panel shows the spectrum obtained when the Csl calorimeter 
veto cut is applied. 

Fig. 14. The predicted PIBETA calorimeter spectra for monoenergetic 69.8 MeV 
positrons and 70.8 MeV photons. ADC values for the Csl crystal with the maximum 
ADC value together with ADC readings of its nearest neighbors were summed. A 
GEANT simulation used different numbers of photoelectrons/MeV and linear axial 
light collection nonuniformities for each crystal, as measured in the RASTA appa- 
ratus. The dashed histograms represent the simulation with the software detector 
gains equal to 1. 

Fig. 15. The RASTA computer screen display showing the fitted histograms during 
automated data acquisition and analysis of a single Csl crystal. 
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Table 1 

Average scintillation properties of the painted hexagonal and pentagonal PIBETA 
Csl calorimeter crystals PENTAs, HEX-As, HEX-Bs, HEX-Cs and HEX-Ds (200 
crystals). The light yields are normalized to the temperature of 18°C. All other 
parameters were measured at the average laboratory room temperature of 22° C. 

BICRON Csl Harkov Csl 
Crystals (22) Crystals (178) 



Fast-to-Total Ratio (100 ns/1 fis gate) 0.835 0.760 

No. photoelectrons/MeV (100 ns ADC gate) 74.2 65.2 

No. photoelectrons/MeV (1 (is ADC gate) 88.8 85.8 

Fast Light Temp. Coefficient (%/°C) -1.20 -1.61 

Total Light Temp. Coefficient (%/°C) -1.35 -1.40 
Axial Nonuniformity Coefficient (%/cm), 

z<10 cm, 100 ns ADC gate 0.178 0.297 
Axial Nonuniformity Coefficient (%/cm), 

z>10 cm, 100 ns ADC gate 0.405 0.635 
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Table 2 

Average scintillation properties of the lacquered half-hexagonal and trapezial PI- 
BETA Csl calorimeter modules HEX-Dl/2s and VETO-l/2s (40 crystals). All 
crystals were polished and painted and wrapped in an aluminized Mylar foil after- 
wards, 

BICRON Csl Harkov Csl 
Crystals (3) Crystals (37) 



Fast-to-Total Ratio (100 ns/1 fis gate) 0.806 0.757 

No. photoelectrons/MeV (100 ns ADC gate) 31.3 26.2 

No. photoelectrons/MeV (1 ^s ADC gate) 38.8 34.6 

Fast Light Temp. Coefficient (%/°C) -2.08 -0.63 

Total Light Temp. Coefficient (%/°C) -1.78 -0.69 
Axial Nonuniformity Coefficient (%/cm), 

z<10 cm, 100 ns ADC gate 0.062 0.315 
Axial Nonuniformity Coefficient (%/cm), 

z>W cm, 100 ns ADC gate 1.592 1.735 
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Table 3 

The predicted energy resolutions and tail contributions for 69.8 MeV e + and 70.8 
MeV 7's events in the full PIBETA calorimeter. The light output in photoelec- 
trons/MeV and the linear axial light collection nonuniformities measured for in- 
dividual Csl crystals (a 2 i i , Z 2 = RASTA) were used in the GEANT simulation. Values 
corresponding to perfect optically homogeneous crystals (a z i tZ 2=0) are shown for 
comparison. 



Parameter 


69.8 MeV e+ 




70.8 MeV 7 




No E v cut 


a>zi,z2=0 


RASTA a zl , z2 


a z i, Z 2=0 


RASTA a zl>z2 


Peak Position (MeV) 


68.87±0.03 


68.78±0.03 


70.10±0.03 


69.50±0.05 


FWHM (2 2o) (MeV) 


3.66±0.03 


3.50±0.03 


3.98±0.03 


3.97±0.04 


5<Events<54 MeV (%) 


6.44±0.09 


6.11±0.09 


8.20±0.10 


7.32±0.09 


5<Events<55 MeV (%) 


6.89±0.09 


6.48±0.09 


8.86±0.10 


7.83±0.10 


E v <5 MeV 




RASTA a z i, z2 


a>zi,z2=0 


RASTA a zl>z2 


Peak Position (MeV) 


68.88±0.03 


68.89±0.03 


70.07±0.03 


69.50±0.04 


FWHM (220 ) (MeV) 


3.67±0.03 


3.49±0.03 


3.99±0.03 


3.97±0.04 


5<Events<54 MeV (%) 


1.91±0.05 


1.70±0.05 


3.92±0.07 


3.15±0.07 


5<Events<55 MeV (%) 


2.24±0.05 


1.99±0.05 


4.46±0.07 


3.57±0.07 
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Table 4 

The predicted energy resolutions and tail contributions for 69.8 MeV e + and 70.8 
MeV 7 events in PIBETA clusters containing the crystal with maximum energy 



CLcpUolLlUIl dllLL Lib IlcdlcoL 


IieigllUUI h. 








Parameter 


69.8 MeV e+ 




70.8 MeV 7 




No E v cut 




RASTA a z i )22 


Q-2i,^2 = 


RASTA a zl:Z2 


Ppak Position fA/TpV^ 


cv-7 74-i-n nQ 

U 1 . 1 'i—L U.Uu 


«7 sn+O 04 


fi8 80+0 03 


so «q+o 03 

UO.UO_LU.UO 


FWHM/™> (MpV\ 

J. VV j-±lvi(]\jp^) I 1V±C V 1 




c Q7+n 04 


« 01 +0 03 


r. qo+o 04 


5<Evpnts<54 MpV (%\ 


U.UOXU.U ( 


q «7+0 07 


c 47+0 08 


4 fi7+0 08 


o<xijvents<oo MeV (Vo) 


/i £ A _i_n n^7 
4.54±U.U / 


A a A _i_n no 
4.o4±U.08 


6.OO±0.09 


o.39±0.09 


E v <5 MeV 


«2l,z2=0 


RASTA a z i ;22 




RASTA a zl>z2 


Peak Position (MeV) 


67.75±0.03 


67.76±0.03 


68.80±0.04 


68.60±0.04 


FWHM (NN) (MeV) 


5.52±0.03 


5.24±0.03 


5.92±0.03 


5.85±0.03 


5<Events<54 MeV (%) 


3.53±0.07 


3.17±0.07 


5.28±0.08 


4.29±0.08 


5<Events<55 MeV (%) 


4.19±0.07 


3.76±0.07 


6.10±0.09 


4.98±0.08 
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